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ABSTRACT 

I present a study of two Herbig Ae stars that are in completely different evolutionary stages: V892 
Tau and PV Cep. Using sub arc-second interferometric observations obtained with the Combined 
Array for Research in Millimeter-wave Astronomy (CARMA) at A=1.3 and 2.7 mm, I have for the 
first time resolved their disks. I deduce that the 5 Myr old V892 Tau has a low dust opacity index 
(9=1.1 and a disk mass of ^0.03M Q . These values correspond to the growth of its dust into large up 
to centimeters size structures. In contrast, the very young (a few xlO yrs) PV Cep has a quite high 
opacity index (9=1.75 and a more massive disk 0.8M Q . PV Cep has the youngest resolved disk around 
any Herbig Ae star. Unlike the youngest T Tauri and Class stars, which contain large and processed 
grains, the young Herbig Ae star, PV Cep, disk contains ISM-like unprocessed dust. This suggests 
that PV Cep's dust evolution is slower than T Tauri stars'. I also present high spatial resolution 
interferometric observations of the PV Cep outflow. The outflow inclination is consistent with the 
orientation of the known Herbig-Haro flow in that region, HH315. 

Subject headings: Stars: formation, pre-main-sequence, circumstellar matter, Herbig Ae/Be, T Tauri 
- Techniques: interferometric 



1. INTRODUCTION 

Grain growth is expected to occur during the pre- 
main-sequence phase of star formation in the circum- 
stellar disks, when the disk is still gas-rich. The dust 
is processed from ISM-like sub-micron size grains into 
mm and cm-size struc tures on their way to form planets 
(jBeckwith et al.l 120001 ) . Interferometric sub-millimetric 
and millimetric observations are playing a key role in the 
detection and study of these disks. 

The intermediate-mass pre-main-sequence, Herbig 
AeBe (HAEBE), stars have a mass ranging roughly be- 
tween 2 to 8M q , with a spectral type range BO to F5 
(|Herbiglll960ft . They are classified in two sub-types. Her- 
big Ae stars that have masses < 5M Q and their spec- 
tral type vary between B5 and F5. Herbig Be stars 
are more massive > 5M Q and earlier type stars, BO- 
BS. They are considered as the progenitors of Vega like 
stars th at are surrounded by debri s disks with large solid 
bodi es (|Backman fc Parescelll993ft and possibly planets 
(e.g. iWilner et al.l 120021 ) . HAEBE stars are still not as 
well studied and understood as their lower mass counter- 
parts T Tauri stars (TTS). Theoretical and observational 
studies of low-mass pre-main-sequence stars (TTS) have 
shown that they form via accretio n from a circumstellar 
disk, with an associated outflow dBeckwith et alJ 11990b 
Mundv et al.lll99"6tlRodmann et al.l l2006: An drews et alJ 
2009). 

From spectral energy observations, iHillenbrand et ail 
(1992) have shown that the infrared excess observed to- 
ward a large sample of HAEBE stars is due to thermal 
emission from the circumstellar material. Using sub- 
millimeter single dish obse rvations on the JCMT with 
the SCUBA array camera, ISandell fc Weintraubl (j2002T ) 
have resolved extended circumstellar emission around 
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two dozen HAEBE stars, and thus suggested the exis- 
tence of massive disks around these stars. High spatial 
interferometric radio observations later played an impor- 
tant role by r esolving circumstellar disks around Herbig 
Ae stars (e . g.lTesti et al.l l2003b lHamidouche et al.l 120061 
llsella et al.l 12007ft . These observations have shown that 
those disks have already processed their dust into larger 
structures. However, the number of known disks around 
Herbig Ae stars is still very small to make absolute con- 
clusions on their evolution. I present a study of two more 
Herbig Ae stars V892 Tau (Elias 3-1) and PV Cep, which 
are in completely different evolutionary status. 

The older sou rce V 8 92 T au was detected first in 
the near-IR by lEliasI (|1978D . V892 Tau is located 
in the Taur us molecular cloud close to the dark 
cloud L1495 (|Skinner et al.l 11993ft. It is of spectral 
type B8±3 ^Hernandez et al.l 120041 ) and of a mass of 
2.8M^ dSmith et all 12005ft at a distance of 130 pc 
(|Torres et al.1 2009f). From its relatively high X-ray 
energ y (lHamidouche et ail 120081 ) and low Ha emis- 
sion (Ma noi et al.l 120061 ) . V892 Tau must be evolved. 



Its age is estimated at 5 Myr (iManoi et al.l _ 
iHernandez et al.l l{20 04) classified it as a main-sequence 
star and hence without a circumstellar dis k. On the other 
hand , from mid-infrared observations, (Monnicr et al. 
|2008|) have resolved it as a close binary (7.7 AU) and 
have shown that it is surrounded by a circumbinary disk. 

The youngest object of this paper is PV Cep. 
iNatta et al.1 (12000ft re ported its age of 0.3 Myr. Using 
the iSiess et alJ (120001) m odel. I deduce an age < lMyr 
(|Hamidouche et al.ll2009l ). It is an embedded star located 
at the head of the cometary nebula GM 29. It is also lo - 
cated slightly below the birthline (|Fuente et al.l fl998) . 
PV Cep has a mass o f 3.5M Q and is classified as an A5e 
star (|The et al.lll994[). In their e v olutio nary classification 
of HAEBE stars. iFuente et al.1 (|1998l ) ranked it in the 
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least evolved type stars, Type 1. PV Cep shows strong 
P Cygni profiles in Balmer li nes ([Hernandez et al.|[2004f) 
and strong and variable Ha (jCohen et al.lll981[) . That is 
due to a strong wind or an outflow. This is consistent 
with the youth of this object and the energetic bipo- 
lar outflow that was dete cted toward this star, extended 
over 2.6 pc in projection ()Levreault| [l984; Reipurt h et al.l 
|1997|) . Spectral energy distribution of P V C ep between 
3-100 /Ltm defines it as a Class I young star (|Stapelfeldtl 
[19961) . 

So far no interfcrometric observations have been done 
toward V892 Tau and PV Cep at 1.3 mm. I present 
the highest resolution aperture synthesis images of V892 
Tau and PV Cep ever done, as well as the first observa- 
tions of PV Cep at 7 mm. In this study, I report that 
V892 Tau and PV Cep are surrounded by circumstellar 
disks and discuss their different disk properties. I focus 
on their dust properties to probe the disk evolutionary 
status by modeling the disk spectral energy distribution 
(SED). Their distinctive disk properties are consistent 
with their respective ages. The dust evolutionary status 
of the youngest PV Cep is different than the youngest 
T Tauri and Class stars. I also present a new inter- 
fcrometric observational evidence of the bipolar outflow 
driven by PV Cep. 

2. OBSERVATIONS 

I made interfcrometric millimeter observations with 
the Combined Array for Research in Millimeter- wave As- 
tronomy (CARMAQ) toward these two Herbig Ae stars 
(Table 1). During my observations period, CARMA was 
a combined array of 15 antennas, six 10.4 m antennas 
and nine 6.1 m antennas. 

I observed V892 Tau at A=1.3 mm in C-array con- 
figuration, which corresponds to a maximum baseline of 
^300 m. These observations were taken in observing 
cycles of 20 min, with 17 min on source and 3 min on 
the phase calibrator. The flux calibrator MWC349 was 
observed at the beginning of the track. In B-array config- 
uration, I observed the star at A=2.7 mm, with a maxi- 
mum baseline of 810 m. The observations were taken in 7 
min cycles, with 5 min on source and 2 min on the phase 
calibrator. The flux calibrator Uranus was observed at 
the beginning of the track. 

PV Cep was also observed at both A=1.3 mm and 2.7 
mm wavelengths with CARMA. In C-array configura- 
tion, I observed PV Cep at A=1.3 mm, with a maximum 
baseline of ~ 300 m. I used observing cycles of 20 min, 
with 18 min on source and 2 min on the phase calibrator. 
The flux calibrator Uranus was observed at the beginning 
of the track. For A=2.7 mm, I used both D- and B-array 
configurations corresponding to maximum baselines of 
110 m and 790 m respectively. I observed in D-array in 
22 min cycles, 17 min on source and 5 min on the phase 
calibrator, and in 8 min cycles, 5 min on source and 3 min 
on the phase calibrator in B-array. For both tracks, the 
flux calibrator MWC349 was observed at the beginning. 

For both sources, I configured the correlator in 
CARMA to have two bands set at 500 MHz for the con- 
tinuum and one band at 8 MHz centered at 230.538 GHz, 
A=1.3 mm, and 115.271 GHz, A=2.7 mm. These config- 
urations provide a total of 1GHz bandwidth in the con- 

1 CARMA has since been combined with the SZA array 



tinuum. I have FFTed and clea ned the data usin g the 
MIRIAD astronomical package (jSault et alj [l995). For 
PV Cep mapping at 2.7 mm, I combined B- and D-array 
observations using robust=-0.5. D-array observations are 
sensitive to the extended circumstellar structures larger 
than 6". The resulting beam size is 0.88"x0.73" with 
a position angle of 22°. In this study, I mainly focus on 
the dust continuum emission. 

I used the VLA (Very Large Array) of the NRAC0 in 
C-array configuration to observe PV Cep with Q-band 
receivers (A=7mm, 43 GHz). The source was not de- 
tected with an rms noise cr=0.19mJy/beam. An upper 
limit to the flux density of 4cr gives 0.76 mJy. 

3. OBSERVATIONAL RESULTS 

V892 Tau - I present in Table 2 the deduced disk pa- 
rameters obtained by fitting a two-dimensional elliptical 
Gaussian from MIRIAD. I estimate the outer radius from 
the Gaussian R g = 1.4 FWHM along the major axis of 
the disk to include 99.9% of the flux, assuming a circular 
shape for the circumstellar disk. I have resolved the disk 
of V892 Tau at A=2.7 mm (Figure 1). Figure 2 shows 
the vector averaged visibility amplitude data binned in 
annuli around the source. The bin size was chosen to 
provide approximately a same number of u, v visibility 
points per bin. The visibility amplitude is plotted as a 
function of the projected baseline length. Compared to 
the point source, I see that the disk is resolved at the 3 a 
level of the amplitude, within error bars. The total flux 
obtained at A = 2.7 mm is consistent wit h low-resolution 
observations of Idi Francesco et al.1 (|1997l ). 

1 deduce a disk inclination of 59° and a P.A. = 64°. 
From mid-infrared imaging of V892 Tau, JM 

onnier et al.l 

2008) deduce a similar value 60° of the disk inclination. I 
deduce a disk radius of R ? ~ 100 AU using a distance of 
130 pc (jTorres et alj|2009l ). At A=1.3 mm, I find a point 
source (Figure 1) with a disk radius smaller than about 
170 AU. The total flux is quite similar wi thin error bars 
with single dish observations of V892 Tau (Henning et al. 
1998). This shows that there is a negligeable contribution 
of an envelope emission in our observations. The flux at 
these wavelengths has a very low contamination by free- 
free emission from the ionized gas 10% and 2%, at 2.7 
and 1.3 mm respectively. 

PV Cep - These CARMA interferometric observations 
reveal circumstellar disk emission around PV Cep at both 
observed wavelengths (Figure 3). The deduced disk pa- 
rameters are given in Table 2. The peak position fits 
the optical position of PV Cep (jDucourant et al.1 120051 ) 
within 0.1", which is ^5% of the beam size. This is 
consistent with the compact rad io emission at 3.6 cm 
detected by lAnglada et al.l (| 19921 ). At A=1.3 mm, I re- 
solved the disk (Figure 3). The total flux of 280 mJy is 
25% higher than the peak. It is lower than single dish 
observations of 360 mJyEI by only 20%. Unlike interfer- 
ometers that filter the surrounding emission, single dish 
observations include emission from the surrounding ma- 
terial. The low difference in flux with single dish observa- 
tions, which have much larger beams > 6" indicates that 
there is an insignificant envelope emission contribution 

2 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation. 

3 G. Sandell private communication 
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to my interferometric observations. Figure 4 is analo- 
gous to Figure 2. It shows that PV Cep is well resolved 
at A = 1.3 mm compared to a point source emission. 

At 2.7 mm, the total flux 3 2 mJy is consisten t with 
low resolution observations of [Stapelfcldt] (|1996| ). Fig- 
ure 4 shows a scattered emission at the longest baselines, 
^250kA. I therefore consider the disk only partially re- 
solved. The disk sizes at 1.3 and 2.7 mm wavelengths 
are consistent. The inclination and the position angle 
are rather different, but the values at 2.7 mm are uncer- 
tain. For the disk modeling in Section 4, I use the values 
from the resolved disk image at 1.3 mm. 

By combi ning existing observa tions of PV Cep at A = 
3.6 cm from lAnglada et al. (1992) with the upper limit at 
A=7mm in this work, I constrain the non-thermal emis- 
sion spectral index a < 0.73 (where S v oc v a ). This 
value is consistent with the value deduced bv lStapelfeldtl 
(1996) ~0.6. A positive value of a means that free-free 
emission is from the ionized gas in the outflow. I deduce 
a negligible free-free contribution to the 2.7 and 1.3 mm 
flux, less than 1.4 mJy and 2.3 mJy respectively. 

From the molecular line CO(l— >0) emission, I resolve 
the inner part of an outflow powered by PV Cep, which 
coincides with the source position and its disk. The out- 
flow position angle is consistent with the disk inclina- 
tion (Figure 5). The NW emission corresponds to the 
blue-shifted part of the outflow at a velocity range of - 
2.3-0.9 km sec -1 . The SE emission corresponds to the 
red-shifted part at higher velocities within the range 4.1- 
8.5 km sec -1 . In my interferometric observations, I only 
see the closest region (~ 0.25 pc) around PV Cep, thus, 
the beginning of the outflow. The corresponding CO 
spectrum is shown in Figure 6. The emission drops at 
2-3 km sec -1 . This corresponds to the emission from 
the cloud, which I resolved out. This valu e is also con- 
sistent with the cloud velocity deduced by ICohen et al.l 
(p8lfl . Fi gure 5 shows a S-shape like outflow emission. 
This can be explained by a rotation of the outflow. This 
suggests that PV Cep outflow rotates around its parallel 
axis counter-clock wise, if seen pole-on from its northern 
side. 

The flow wide opening angles show that the out- 
flow swept the sur r ounding material of t he cloud 
(jReipurth et alJU997t lArce fc Goodmar3l2002l ). This is 
likely why I clearly detect the circumstellar disk around 
the very youn g embedded PV Cep star. Using single dish 
observations, lArce fc Goodm an (2002) mapped a large 
region (a few pc) surrounding PV Cep. They find a not 
well collimated outflow in the north of PV Cep. This 
is consistent with the wide angle of the V-shaped blue 
shifted part (Figure 5). 

4. DUST PROPERTIES 

The interstellar dust is processed in the flared disk or 
atmosphere component. The grown and large grains sed- 
iment toward the disk midplane immediately while the 
interstellar micron size d ust may take abou t a Myr to 
reach the midplane (e.g. iNatta et al.l |2000| ). The disk 
midplane contains the bulk of the material and this is 
where the grains grow into larger structures and plan- 
etesimals by collisions. 

4.1. Disk Parameters and Model 



For this study, I model the circumstellar disk midplane 
as an optically thick disk and geometrically thin with 
a hole in the middle. The disk midplane emission 
is dominant from the far-infrared into the millimeter 
region of the SED, or the Rayleigh- Jeans component. 
Modeling this component constrains the dust prop- 
erties relevant for this study. I add 20% calibration 
uncertainties to the systematic errors for my new 2.7 
and 1.3 mm CARMA data. In the sub-millimeter, 
I use new data from SCUBA at 450 and 850 /im. 
(ISandell. C. and Weintraub. D. A., and Hamidouche. MJ 
20101 submitted to ApJS.). At A=7 mm, I use PV Cep 
flux as an u pper limit, and for V892 Tau, I use the flux 
value from Idi Francesco et all (|1997f ). I also add the 
existing MIPS3 A=70/zm observations of V892 Tau. 

I model the disk midplane withou t an envelope and 
using a power-law chara c terization dMundv et al.l 119961 : 
lHamidouche et all [20061 : lAndrews fc Williams! I2007alf . 
This is a simple and convenient model for this study. 
The disk midplane emission is presented by 

fRoui —Y,(R)k R 

F„ = 27roost / {l-exp[ y -^}}xB v [T d {R)] —dR 

J R . cos i 

(1) 

The dust opacity is characterized by the power-law 
k(v) = Ki2oo(i / /1200) _/3 , where v is the frequency, j3 is 
the opacity index, and K1200 — 0.01 cm 2 g _1 . This in- 
cludes the gas contribution to the t otal disk mass wit h 
a gas-to-dust mass ratio of 0.01 (e.g. INatta et alJ feOOO). 
Due to the low-spatial resolution of my observations, I 
do not use a radial profile for the opacity index (3. The 
disk mass Mp is deduced from the constrained value of 
Mi?xki2oo- The distance to the star is D. The disk ra- 
dius R varies from the disk inner radius, which is the 
middle hole radius, Ri n to the disk outer radius R ou t- 
The surface density profile is defined as = Y^iau x 

R~ p . For the temperature profile in the disk midplane, I 
use a radial power-law T(R) = T\au x R~ q characteri- 
zation, which is estimated by lAndrews et al.l ()2009l ) . 

The disk outer radius and surface density distribu- 
tion cannot be constrained from the SED model ing (e.g . 
lHamidouche et all [20061 : lAndrews fc Williams! I2007bf t. 
The spatial resolution provided from these observations 
is not sufficient (Figure 2 & 4) to constrain the sur- 
face density profile and therefore the outer radius from 
the disk image and visibility. I fix the outer radius 
Rm/,t = Rn, whi ch is the best estimate here (see also 
iRicci et"aTI l2010l ). I fix the value of p=l that is the 
value found i n the well resolved disks around Herbig 
Ae stars (e.g.lTesti et al.l l2003t lHamidouche et al.l 12006 1 : 
llsella et al.ll2007| ). Nevertheless. I tested the SED model 
by varying R ou t and for p=0.5-1.5, where 1.5 value is the 
derived value for the s olar-nebula surface density distri- 
bution (jHavashil 119811 ). I found that p values do not 
significantly affect the two relevant disk parameters val- 
ues, within error bars, the opacity index j3 and the disk 
mass Me>. 

The inner rad ius is not relevant for this study 
(|Testi et al.l 120031 ). So I fix its value to the sublimatio n 
radius of Herbig Ae stars i? m =0.3 AU (|Vinkovicll20"06T ). 
The disk inclination is obtained in the previous section. 

4 Spitzer Campaign Id 749 (MIPS004000) 
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I keep Tiau, the index q, the dust opacity index /3, and 
disk mass Mjj as free parameters. I deduce the models 
that fit best the observations using a x 2 minimization 
method for the first three parameters spanning a three- 
dimension cube. Once the best fit 3D grid is found, the 
disk mass is adjusted until I find the lowest x 2 ■ I eval- 
uate the models that are close to the minimum x 2 by 
taking the models with a likelihood > 0.5. I discuss in 
the following section the dust opacity and disk mass. 

4.2. Modeling Results and Discussion 

I made well sampled SEDs from the far-infrared to A=7 
mm. The long wavelengths allow probing the dust prop- 
erties to centimeter size. I deduced from the disk model 
the dust opacity index f3. I calculated the millimeter 
spectral index a mm for both stars at millimeter wave- 
lengths. I used only the interferometric data (1.3, 2.7 
and 7 mm) to avoid any over estimation of the surround- 
ing cloud emission, and thus the calculated a mm value is 
robust. 

When the emission is optically thin, the spectral in- 
dex is well correlated with the opacity index a rnm -2=/3 
(jBeckwith fe Sargent] Il99ll ). An optically thin emission 
probes directly the dust properties independent of a ny 
complexity of the disk morphology (iTes ti et aLll2003ft . I 
varied the value of /? in my model from to 2. They 
correspond to a disk cont aining only non-processed in- 
terstellar like dust, f3 >1.7 (Wcingart ner fc DrainduOOlft 
or dust comp letely grown into very large particles ( "peb- 
bles"), /3=0 (jBeckwith fc Sargentill99ll ). An intermedi- 
ate value j3 ~1 is consistent with t he strong contribu - 
tion of the gr ains larger than 1 mm (iNatta et al.| [2007). 
Furthermore, iMivake fc Nakagawa ( 1993) interpret such 
values ~1 corresponding to dust sizes from 0.1-10 cm. 
The results of SED modeling are presented in Table 3. 

V892 Tau - SED is shown in Figure 7. A typical fit 
model is shown in a solid line. Deprojected visibility 
profile of this typical fit model in Figure 2 (solid line) 
demonstrates qualitatively that the model and visibility 
data are consistent to within errors. I deduce from the 
best fit models an opacity index /3= 1.1 ±0.1 and a disk 
mass of 0.035±0.015 M Q (Table 3). The corresponding 
disk-to-star mass ratio Md/M* = 0.01, or a disk mass 
of 1% of the stellar mass. The low value of /3 shows that 
the disk has passed its initial stage and processed the in- 
terstellar dust. This is consistent with the disk low mass, 
since the large structures' emission cannot be seen due to 
their small surface area. I do not find any model that fits 
the observations with an opacity /3—0. Thus, the dust is 
not completely processed into "pebbles". I deduce the 
millimeter spectral index a mm =2.93, and thus the ra- 
tio of optically thick to optically thin emission is only 
a few percent at millimeter wavelengths. The contribu- 
tion of the inner optically thick region is negligible. This 
shows that my results are quite robust and the shallow 
observed slope in the SED is due to processed grains. 
The disk morphology does n ot affect the deduced grain 
properties (|Testi et all 120031 ). 

V892 Tau is a relatively old pre-main-sequence star sur- 
rounded by an evolved disk. V892 Tau has gone through 
a large step of evolution. Furthermore, the V892 Tau 
disk's evolutionary phase is similar to the Herbig Ae 
star's, MW C 275 (HP 163296). T hey have the same 
age, 5 Myr (jMontesinos et al.l 120091 ). and a similar mass 



~2.5M (jlsella et al.ll2007D . They show a shallow mil- 
limeter slope with a mm ~3 and an optically thin emis- 
sion with an opacity index /?=1. Nevertheless, MWC 27 5 
has a different spectral type, A3 (jBlondel fc Diid l2006'). 
and is isolated. In contrast, V892 Tau is an earlier type 
star (B8) and located in Taurus. The similarity of their 
disks' evolutionary stage shows that their evolution is not 
correlated with their spectral type. They also have simi- 
larHa EW ~20 A, and thus accretion rate (jManoi et al.l 
120061 and references there in) . Their environment clearly 
does not affect their disk evolution or accretion process. 

PV Cep - The SED is shown in Figure 8. The solid line 
shows a typical fit disk model. As in Figure 2, the corre- 
sponding deprojected visibility profile in Figure 4 (solid 
line) demonstrates qualitatively that the model and visi- 
bility data are consistent to within errors. I deduce from 
the best fit models 0=1.75 and M D =0.76M Q (Table 3). 
The corresponding disk-to-star mass ratio Md/M, =0.2, 
or a disk mass of 20% of the stellar mass, which is a rela- 
tively high value compared to older Herbig Ae stars. The 
dust is clearly not processed into large structures and is 
rather similar to the ISM (sub-)micron size dust. It has 
a mm spectral index a mm =3.58. I deduce an optically 
thick emission contribution of only few percent at mil- 
limeter wavelengths. The youth of the star < lMyr is 
consistent with the disk evolutionary stage. In fact, PV 
Cep was classified as an A5 star, but, equivalent to a 
very young T Tauri star that is c oming out of its prestel- 
lar envelope ()Cohen et alJ ll981). This is the youngest 
resolved disk around any Herbig Ae star. 
The star has reached the pre-main-sequence stage with- 
out processing its dust. In contrast, the lower mass T 
Tauri stars appear to have already processed their dust 
at their youngest age. Compared to PV Cep, interfero- 
metric observations of eight young (a few x 10 5 yrs old) 
TTS provided mu ch smaller dust opacity indexes /3 <1 
(jlsella et al.1 12009ft. From a stud y of two dozen TTS of 
different ages. IRicci et al.l (|2010ft have found that all the 
youngest stars < 1 Myr have j3 < 1. Similarly, interfero- 
metric observations of three Class sources have shown 
that t heir grains are alre ady processed and grown at this 
stage (|Kwon et al.l l2009) . Additionally, even the massive 
(IOMq) main-sequence star MWC 297 has a very low- 
mass disk with evolved grains /3=0. 1-0.3 which is pre- 
sumably remnant from a more mas sive disk when the 
star was younger (|Manoi et~aT][2007l) . If the Herbig Ae 
star, PV Cep, would evolve to a main-sequence star in 
a similar way as TTS, it should have already processed 
its dust into large grains. This makes the PV Cep case 
important and worth further investigations. 

I suggest that the disk of PV Cep will evolve and the 
dust grow into large structures. Since PV Cep is rela- 
tively massive M*=3.5M© compared to TTS, its hor- 
izontal gravity force, GM^m,d/r 2 , acting on the orbit- 
ing dust, rrid, at a distance r in the disk atmosphere is 
stronger. The strong stellar horizontal gravity may hold 
the grains for a longer time in the PV Cep disk atmo- 
sphere th an in TTS's, before th ey sediment toward the 
midplanc (Be ckwith et a l. 2000). Nevertheless, while ro- 
tating i n the disk atmosphere, the dust collides and coag- 
ulates (|Dullemond et al.l 12007ft . and thus becomes mas- 
sive enough to be dragged down toward the midplane by 
the combined disk and stellar vertical forces (|Natta et al.l 
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120071) . At that time, the dust would grow further in 
the midplane via collisions or local gravitational ins tabil- 
ity (|Goldreich fc Wardl[l97l lYoudin fc ShTJ l200l . and 
thus show similar dust properties as in the more evolved 
Herbig Ae disks, such a s V892 Tau's and MWC 480's 
(jHamidouche et all 120061 ). Additionally, the longer du- 
ration of the dust evolutionary process in the young PV 
Cep disk compared to the youngest T Tauri stars' ap- 
pears to be independent of the disk mass or its gravi- 
tational force. In fact, PV Cep disk-to-star mass ratio 
(20%) is similar to the average T Tauri stars' ratio. 

5. CONCLUSION 

I have reported in this paper a study of two 
intermediate-mass Herbig Ae stars that are very different 
in their evolutionary stages and were born in two differ- 
ent clouds. I studied the early type B8 and quite evolved 
5 Myr old star V892 Tau and the later type A5 and 
young a few xlO 5 yr deeply embedded PV Cep. I used 
new CARMA interferometric observations at 1.3 and 2.7 
mm to resolve their disks and look into their dust. I have 
shown that their dust properties are completely different 
and this is consistent with the two different stellar evo- 
lutionary stages. The main results of this study are: 

1) I made the highest resolution interferometric mil- 
limeter observations ever made of the disks of V892 Tau 
and PV Cep, resolving for the first time their disks. They 
are both optically thin at A = 1.3 mm. 

2) I deduce that V892 Tau is an evolved pre-main- 
sequence star. It has passed through an important phase 
of its evolution and has a dust opacity index 0=1.1 and 
disk mass 0.035M©, which is about 1% of the stellar 
mass. The dust in the disk are grown into large cen- 
timeter size structures. I find an interesting resemblance 
between V892 Tau and another 5 Myr old isolated Her- 
big Ae star MWC 275. From their low disk mass and 
opacity index (5 ~ 1, they seem to evolve in a similar 
way independent of their spectral types and surrounding 
environment. 

3) I show interferometric observations of the CO(l— ?-0) 
molecular line emission toward PV Cep. I see clearly an 
outflow driven by PV Cep that is consistent with its disk 
inclination. The outflow inclination is also consistent 
with the orientation of the giant known Herbig-Haro HH 
315 flow in that region. The observations suggest that 
the outflow is rotating around its parallel axis. 
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4) PV Cep is a very young pre-main-sequence star that 
just left its prestellar envelope. It is surrounded by a 
massive disk of 0.76M Q , which is relatively massive for 
a Herbig Ae stars. It is 20% of the stellar mass. The disk 
is composed of unprocessed ISM-like dust, corresponding 
to the opacity index /3=1.75. Unlike Class or young 
T Tauri disks, the PV Cep disk does not contain large 
structures or grown grains. I suggest that due to the 
higher mass of PV Cep compared to T Tauri stars', the 
grains in the PV Cep disk atmosphere will take a longer 
timescale to grow and settle downward toward the disk 
midplane. At that time, they will presumably continue 
their evolution like the more evolved Herbig Ae stars, 
such as V892 Tau. The PV Cep star is the youngest and 
the least evolved resolved disk known around any Herbig 
Ae star. Its dust evolutionary stage is different than the 
youngest T Tauri and Class stars'. 

This raises the question of whether the evolutionary 
process of Herbig Ae stars is different than T Tauri stars'. 
In this study, PV Cep case indicates that it is different, 
or slower, in the beginning, at the age of a few x 10 5 yrs. 
Although, a study of more disks around more Herbig 
stars of different ages is important to further probe 
the evolution of Herbig Ae compared to T Tauri stars' 
(lHamidouche. M.. and Loonev. L. W.. and Mundv. L. G.I 
120101 in preparation.). Higher spatial-resolution in- 
terferometric observations are necessary for such a 
study. The new generation of interferometers ALMA 
and EVLA will provide sensitive high-spatial resolution 
images. This will allow us to extend the sample of 
resolved Herbig AeBe disks into the fainter and farther 
objects, such as the massive Herbig Be stars. 
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Fig. 1. — CARMA images of V892 Tau. The left image is at A=2.7 mm and the right at A=1.3 mm. The contour steps and the starting 
levels are shown in each image. 
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Fig. 2. — V892 Tau - Vector averaged visibility amplitude as a function of projected baseline length at A=2.7 mm. The error bars are 
standard deviation of visibilities in each bin. This u, v amplitude decreases with the baseline length down to 3 a amplitude level within error 
bars. This suggests that the disk is resolved. The dashed line corresponds to a point source emission model and is shown for comparison 
with the disk emission. The solid line (with open squares) shows deprojected visibility profile from the model (Section 4.2). This is a 
qualitative demonstration that the model and visibility data are consistent to within errors. 
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Fig. 3. — CARMA images of PV Cep. The left image is at A=2.7 mm and the right at A=1.3 mm. The contour steps and the starting 
levels are shown in each image. 
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Fig. 4. — PV Cep - Vector averaged visibility amplitude as a function of projected baseline length. The error bars correspond to the 
standard deviation of visibilities in each bin. (a) A = 2.7 mm - The disk is not a point source. But, because of the decoherence at the 
longest baselines, 250kA, I do not consider the disk fully resolved, (b) A = 1.3 mm - The data amplitude decreases with the baseline length 
down to 3 a amplitude level within error bars. This suggests that the disk is resolved. The dashed lines correspond to a point source 
emission model and are shown for comparison with the disk emission. Similar to Fig. 2, for qualitative demonstration, dcprojcctcd visibility 
profiles arc plotted with solid lines (with open squares) in both figures, showing consistency of model and data visibilities to within the 
errors (Section 4.2). 
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Fig. 5. — PV Cep CO(l— >0) molecular line emission map. The outflow emission is shown in contours for the blue- and red-shifted parts. 
The ellipse marks the disk position. This emission is obtained after subtracting the continuum emission and convolving the image by a 8" 
beam. The contours are in steps of 0.12 and 0.14 K X km sec -1 for the red and the blue respectively. 
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Fig. 6. — The spectrum of PV Cep from CO(l— >0). The channel width is 0.62km sec 1 . The double peaks are consistent with an outflow 
emission. The emission drops at 2-3km sec -1 , which corresponds to the cloud velocity. 




Fig. 7. — SED of V892 Tau after correcting the fluxes from free-free contamination. The solid-line is a typical fit obtained for 0=1, 
M D =0.028M Q , g=0.8, and T 1AU = 550 K the corresponding \l -0-9. CARMA 1.3 and 2.7 mm data are shown with open circles. The 
millimeter spectral index a mm is shown in the top right, and F oc A -2 93 . 
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Fig. 8. — SED of PV Cep after correcting the fluxes from free-free contamination. The solid-line is a typical fit obtained for /3=1.7, 
M£>=O.72M0, <jr=0.5, and Tiau = 700 K the corresponding \r —0.8. The CARMA data at 1.3 and 2.7 mm data are shown with open 
circles. The VLA data at 7 mm is an upper limit. The millimeter spectral index a mm is shown in the top right, and F <x A -3,58 . 
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TABLE 1 

Interferometric Observations. 



Object 


RA 


DEC 


A 


Array 


Obs. 


Beam 


Beam 


Phase 




(J2000) 


(J2000) 


(mm) 


Config. 


Date 


FWHM(") 


P.A.(°) 


Calibrators 


V892 Tau 


04 18 40.60 


+28 19 16.7 


1.3 


CARMA C 


2007 Sept. 07 


1.28 x 0.92 


-85 


3C111 








2.7 


CARMA B 


2008 Jan. 14 


0.87 x 0.53 


-61 


3C111 


PV Cep 


20 45 53.96 


+67 57 38.9 


1.3 


CARMA C 


2007 Sept. 19 


1.07 x 0.74 


-59 


1849+670 






2.7 


CARMA D 


2007 Jun. 06 


6.4 x 4.8 


65 


1927+739 








2.7 


CARMA B 


2008 Feb. 14 


0.85 x 0.67 


22 


1927+739 








7 


VLA C 


2005 Jul. 25 


0.65x0.54 


-8.4 


2006+644 



TABLE 2 

Continuum Observations Results 



Object 


A 

(mm) 


Size 
(") 


incl. P.A. 

(°) (°) 


Rg 

(AU) 


Total Flux 
(mjy) 


Disk 
Resolved: 


V892 Tau 
PV Cep 


1.3 
2.7 
1.3 
2.7 
7 


0.51x0.27 
0.68x0.32 
0.64x0.47 


59+6 64+5 
62+4 -63+3 
43+15 12+14 


< 170 
93 
480 
450 


251 

39 
280 

32 
<0.76 


No 
Yes 
Yes 
Partially 


TABLE 3 

Disk properties deduced from SED modeling 


Object 


M D 

(M e ) 




/3 Tiau 
(K) 


q 


Xr 




V892 Tau 
PV Cep 


0.035+0.015 
0.76+0.04 




1.1+0.1 550+50 
1.75+0.05 700 


0.8 
0.5 


<0.9 
<0.9 





